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ABSTRACT 


This  thesis  presents  an  analysis  of  atmospheric  refrac- 
tivity  conditions  in  the  arctic  Marginal  Ice  Zone.  The 
fundamental  principles  of  atmospheric  effects  on  electromag¬ 
netic  wave  propagation  are  presented  along  with  methods  to 
assess  these  effects.  Using  these  principles,  a  determina¬ 
tion  of  the  occurrence  of  atmospheric  ducts  is  made  from 
meteorological  data  gathered  during  the  Arctic  Marginal  Ice 
Zone  Exercise,  1983  (MIZEX-83) .  The  optimum  coupling  height, 
thickness  and  intensity  for  all  ducts  are  calculated  and 
statistically  analyzed.  Duct  occurrence  is  related  to  area 
meteorological  conditions.  Based  on  area  ref ractivity ,  the 
potential  for  enhanced  area  electronic  warfare  operations  is 


evaluated 
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X .  INTRODUCTION 


A.  BACKGROUND 

Modern  Naval  warfare  will  rely  heavily  on  command  and 
control  communications  (C^),  command  and  control  communica¬ 
tions  countermeasures  (C^CM)  as  well  as  weapons  guidance  and 
other  electronic  warfare  support  systems.  All  systems  em¬ 
ployed  for  these  purposes  propagate  electromagnetic  energy 
in  the  electromagnetic  frequency  spectrum  (Table  1),  and  all 
are  affected,  to  some  extent,  by  the  meteorological  environ¬ 
ment  in  which  they  are  operated.  A  more  detailed  descrip¬ 
tion  of  the  electromagnetic  spectrum  can  be  found  in  Tipler 
(1976).  Task  force  commanders  will  be  able  to  gain  a 
decided  advantage  if  they  can  determine  which  systems  will 
be  enhanced  and  which  will  be  degraded  and  use  this  informa¬ 
tion  in  making  their  tactical  decisions. 

Various  projects  to  study  meteorological  conditions  and 
document  atmospheric  anomalies  that  affect  electromagnetic 
wave  propagation  have  been  conducted  in  the  more  common 
Naval  operating  areas  of  the  world  such  as  the  Mediterranean 
and  Indian  oceans.  Throughout  the  world,  meteorological 
data  are  gathered  daily.  A  comprehensive  statistical  analy¬ 
sis  of  this  data  has  been  conducted  by  Leigh  N.  Orthenburger 
of  GTE  Sylvania  (Orthenburger,  et  al.,  1978).  Meteorologi¬ 
cal  conditions  in  the  area  around  the  northern  ice  pack 
have  just  recently  gained  interest.  An  experiment 


Ill . 


ATMOSPHERIC  EFFECT5  ON  ELECTROMAGNETIC 
WAVE  PROPAGATION 


The  atmosphere  contains  two  specific  regions,  the 
troposphere  and  the  ionosphere,  that  have  a  dramatic  effect 
on  electromagnetic  wave  propagation.  The  ionosphere  extends 
from' approx  imately  31  km  to  400  km  above  the  Earth's  surface 
and  is  composed  of  layers  of  ionized  molecules  that  refract 
electromagnetic  waves  in  the  high  frequency  portion  of  the 
spectrum  (3  to  30  MHz).  Between  the  surface  of  the  Earth 
and  the  upper  atmosphere,  extending  up  to  about  12  km,  is 
the  troposphere.  Refractive  conditions  in  the  troposphere 
primarily  affect  signals  at  frequencies  above  the  high 
frequency  band.  The  troposphere  is  an  inhomogeneous 
dielectric  medium  made  up  of  various  constituents. 
Concentrations  of  permanent  gases,  such  as  oxygen  and 
nitrogen,  remain  relatively  constant;  however, 
concentrations  water  vapor  and  aerosols  such  as  suspended 
organic  and  inorganic  matter  vary  greatly.  Also,  fog,  rain, 
snow,  smoke  and  dust  can  frequently  be  present.  These 
molecules  all  have  their  particular  resonant  frequencies  and 
dipole  characteristics.  Electromagnetic  waves  coming  in 
contact  with  or  in  close  proximity  to  these  molecules  can  be 
absorbed,  scattered  or  refracted  depending  upon  the 
frequency  of  the  electromagnetic  wave  and  the 
characteristics  of  the  molecule.  As  such,  electromagnetic 
waves  transmitted  in  or  entering  the  troposphere  are 


T 


y  *1!  -  \  W-J  -H.  "T 


FIGURE  5 

Continuous  Refraction  of  a  Wave  when  Encounterir 
a  Refractive  Gradient  of  Successive  Dielectric 
Layers  with  Increasing  Indices  of  Refraction 
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the  angle  of  refraction  6  c  equal  to  90  degrees  (Tipler 
197S)  . 


SIN  S c  =  n2/ni  C5.6) 

Thus  far,  only  the  case  of  a  wave  encountering  a 
single  discontinuity  in  the  index  of  refraction  has  been 
considered.  However,  if  successive  dielectric  layers,  each 
constituting  a  discontinuity,  are  stacked  one  upon  the 
other,  a  refractive  gradient  can  be  approximated  as  shown  in 
figure  5.  As  these  layers  become  differentially  small,  the 
wave  can  be  refracted  ,  in  the  example  in  figure  5,  down¬ 
ward.  Further,  if  the  wave  is  propagated  within  the  criti¬ 
cal  acceptance  angle,  it  can  be  trapped  within  the  gradient. 

The  basic  concepts  of  electromagnetic  wave 
propagation  in  a  homogeneous  dielectric  medium  can  also  be 
used  to  describe  propagation  in  an  inhomogeneous  dielectric 
medium  like  the  atmosphere.  However,  unlike  a  homogeneous, 
isotropic  dielectric  medium,  the  atmosphere  can  be 
characterized  by  continuously  changing  values  of  refractive 
gradient  throughout  a  given  vertical  cross  section.  This 
property  of  the  atmosphere  where  refractive  index  changes 
with  altitude  gives  rise  to  atmospheric  phenomena  having 
unique  effects  on  electromagnetic  wave  propagation.  These 
effects  and  the  meteorological  conditions  that  cause  them 
will  be  described  in  chapter  III. 


FIGURE  4 


Total  Internal  Reflection 
Critical  Angle 


toward  the  greater  n. 


Air 


in  this  example,  has  an  index 


of  refraction  which  limits  the  speed  of  the  propagating  ray 
to  a  velocity  vi.  Likewise,  the  water  has  a  different  index 
of  refraction  which  limits  the  speed  of  the  ray  to  a  slower 
velocity  V2-  Consequently,  distances  traveled  by  the  ray  in 
the  two  media  in  the  same  time  are  not -the  same,  resulting 
in  a  change  in  the  direction  of  the  wave  as  the  wavefront 
passes  through  the  discontinuity  boundary.  The  following 
equation,  which  is  Snell's  law,  can  be  derived  by  applying 
geometric  arguments  to  the  drawing  in  figure  3  (Tipler 
1976)  . 


ni  SIN  0 l  =  n2  SIN  0 2  C2.4) 

Equation  <2.12)  can  be  solved  for  0 2  yielding 

SIN  02  =  ni/n2  SIN  d\  (2.5) 

If  n2  is  greater  than  nj,  equation  <2. 5)  can  be  solved  for 
the  angle  of  refraction  02  f°r  sny  angle  of  incidence  0 i . 
If  n2  is  less  than  ni,  equation  <2. 13)  gives  values  for  sine 
greater  than  1  if  the  angle  of  incidence  is  greater  than 
some  critical  angle  0C-  Since  the  value  of  the  sine  func¬ 
tion  cannot  be  greater  than  1,  this  implies  that  there  is 
an  angle  of  incidence,  beyond  which,  the  wave  will  be  total¬ 
ly  reflected  into  the  incident  medium.  The-  critical  angle, 
as  shown  in  figure  4,  is  defined  as  the  angle  which  makes 
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FIGURE  3 


Refraction  of  Plane  Waves  at 
the  Surface  of  an  Index  of 
Refraction  Discontinuity 


Description  of  Refraction 


The  ratio  of  the  speed  of  an  electromagnetic  wave  in 
free  space  to  that  in  a  dielectric  medium  is  known  as  the 
absolute  index  of  refraction  n.  The  index  of  refraction  is 
also  equal  to  the  square  root  of  the  dielectric  constant 
(Hecht  and  Zajac,  1979) . 


n  =  C/V  =  (K)1''2 


(2.3) 


Where: 

C  -  The  speed  of  light 
V  *  Speed  in  the  medium 

Since  the  speed  of  propagation  in  a  dielectric  medium  cannot 
be  greater  than  that  in  free  space,  n  will  always  be  greater 
than  or  equal  to  1.0.  According  to  Hecht  and  Zajac  (1979), 
the  index  of  refraction  for  air  varies  between  1.0003  and 
1 . 0005 . 


When  an  electromagnetic  wave  encounters  a  discon¬ 
tinuity  in  n  like  a  ^ ight  ray  going  from  air  to  water,  the 
path  of  the  incident  wave  will  be  bent  or  refracted.  Figure 
3  depicts  a  light  wave  passing  from  air  into  water.  Since 
light  propagates  as  an  electromagnetic  wave,  a  wave  front  is 
used  here  to  illustrate  the  refractive  process.  Similar 
refractive  conditions  can  exist  between  any  two  dielectric 
media  with  differing  indices  of  refraction,  such  as  two 
differing  air  masses.  As  the  index  of  refraction  for  water 
is  greater  than  for  air  (1.33  as  opposed  to  1.0003),  n2  is 


greater  than  ni  and  it  can  be  seen  that  the  wave  bends 


dipole  moment  as  a  result  of  unequal  sharing  of  valence 
electrons  (Figure  2a) .  When  a  dielectric  medium  is  subject 
to  an  applied  electric  field,  the  internal  charge  distribu¬ 
tion  distorts  under  its  influence.  This  corresponds  to  the 
generation  of  electric  dipole  moments.  These  dipole  moments 
experience  a  torque,  which  tends  to  align  them  in  the  direc¬ 
tion  of  the  field  (Figure  2b) .  This  alignment  of  the 
molecular  dipoles  produces  an  additional  electric  field  due 
to  the  dipoles  which  is  in  the  direction  opposite  the  origi¬ 
nal  field.  The  original  field  is  weakened  in  this  way.  Even 
if  the  molecules  of  the  dielectric  medium  are  nonpolar,  they 
will  experience  an  induced  dipole  moment  (in  the  direction 
of  the  original  field)  in  the  presence  of  the  electric  field 
of  the  propagating  wave.  Here  again,  the  additional  elec¬ 
tric  field  due  to  these  induced  moments  weakens  the  original 
field  (Tipler,  1978). 

The  interaction  of  a  dielectric  medium  with  a 
propagating  electromagnetic  wave  is  described  by  the 
dielectric  constant  K  of  the  medium.  If  Eo  is  the  original 
field  without  the  dielectric,  the  new  field  E  is 

E  =  Eo/K  (2.1) 

The  dielectric  constant  is  the  ratio  of  the  magnitude  of  the 
field  in  free  space  to  the  field  in  the  dielectric  medium. 
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FIGURE  1 

Electric  and  Magnetic  Field 
Vectora  in  a  Plane  Electromagnetic  Wav 
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II.  ELECTROMAGNETIC  WAVE  PROPAGATION 

A.  PROPAGATION  IN  FREE  SPACE 

James  Clerk  Maxwell  demonstrated  that  electromagnetic 
waves.  in  free  space.  propagate  in  a  straight  line  at  the 
speed  of  light  <Xipler,  1976). 

Maxwell  showed  that  electromagnetic  wave  propagation  can 
be  graphically  displayed  on  a  three  dimensional  (XYZ)  plot 
as  shown  in  figure  1.  Note  that  the  electric  field  <Ey>  and 
the  magnetic  field  (Bz)  are  orthogonal  to  each  other  and 
that  the  direction  of  propagation  is  parallel  to  the  vector 
<E  X  B)  and  is  in  the  X  direction. 

This  property  of  propagation  in  a  straight  line  and  at 
the  speed  of  light  is  common  to  electromagnetic  waves  of  all 
wavelengths  as  long  as  the  propagation  path  is  in  a  vacuum. 
Electromagnetic  waves  can  also  be  propagated  in  media  such 
as  water.  glass  and  the  Earth's  atmosphere.  However,  when 
this  happens,  the  propagation  path  can  be  altered  by  the 
medium . 

B.  EFFECTS  OF  A  DIELECTRIC  MEDIUM  ON  PROPAGATION 

1 .  Description  of  a  Dielectric  Medium 

When  an  electromagnetic  wave  is  transmitted  in  a 
homogeneous,  isotropic  dielectric  medium,  the  medium  tends 
to  resist  the  wave  and  the  net  effect  is  to  reduce  the 
propagation  speed.  The  medium  can  be  modeled  as  consisting 


of  numerous  randomly  oriented  molecules  which  have  a 


electromagnetic  wave  propagation  exist  and  to  correlate 
these  findings  to  the  overall  meteorological  conditions  in 
the  area.  This  "area  refractivity  summary"  can  then  be 
used  by  tactical  commanders  operating  in  the  area. 

C.  ORGANIZATION 

The  first  portion  of  this  thesis  (Chapters  II  and  III) 
describes  the  basic  concepts  necessary  to  understand  the 
data  analysis  process.  Chapter  II  deals  with  basic  electro¬ 
magnetic  wave  propagation  in  free  space  and  how  these  waves 
can  be  refracted  in  a  generic,  homogeneous  dielectric  medium 
such  as  glass  or  water.  Chapter  III  covers  the  atmosphere. 
Atmospheric  refractivity  is  introduced  and  related  to  atmos¬ 
pheric  temperature.  pressure  and  humidity.  The  concept  of 
changing  refractivity  with  altitude  or  refractivity  gra¬ 
dients  is  presented  here  and  it  is  shown  how  to  use  these 
gradients  to  detect  atmospheric  anomalies  that  can  affect 
propagation.  Chapter  III  also  describes  refractivity 
assessment  methods. 

The  final  portion  of  this  thesis  deals  with  data  analy¬ 
sis  and  conclusions.  In  chapter  IV,  data  acquisition  and 
analysis  methods  are  described.  A  compilation  of  the  data 
is  presented  along  with  results.  Chapter  V  consists  of  an 
area  refractivity  summary,  an  application  of  this  summary  to 
electronic  warfare  and  recommendations.  Data  are  contained 


in  the  appendices. 


was  begun  in  the  summer  of  1983  to  determine  and  evaluate 
those  meteorological  factors  which  significantly  affect  or 


are  affected  by  the  position  of  the  ice  e> je.  The  Marginal 
Ice  Zone  Experiment,  1983  (MIZEX-83)  is  the  first  large 

A  -■ 

scale  effort  to  assess  meteorological  conditions  in  the  I 

northern  marginal  ice  zone  <MIZ> .  Now  that  interest  and  the 
level  of  research  in  the  MIZ  has  risen,  the  U.S.  Army  Cold 
Regions  Research  and  Engineering  Laboratory  publishes  the  I 

MIZEX  Bulletin  series  which  consists  of  articles  on  the 
results  of  meteorological  and  oceanographic  research  being 

fe  i  > 

done  in  the  area.  Interest  in  the  MIZ  has  led  to  the  | 

development  of  an  area  atmospheric  boundary  layer  model  by 
J.E.  Overland  (Overland,  et  al.,  1983)  and  a  wind  drag 
f|  coefficient  study  by  3. A.  Macklin  (Macklin,  1983).  These  ^ 

and  other  related  reports  can  be  found  in  the  March  30,  1983 
issue  of  the  Journal  of  Geophysical  Research  which  is  a 
H  special  issue  on  the  marginal  ice  zone.  Data  gathered  i 

during  MIZEX-83  will  be  the  basis  of  this  thesis. 

B.  PURPOSE 

I 

Thus  far,  research  on  the  MIZ  has  dealt  with  the  clima¬ 
tology  of  the  area.  Minimal  work  has  been  done  to  correlate 
area  meteorological  conditions  and  atmospheric  ref racti vity . 

» 

The  purpose  of  this  thesis  is  to  analyze  meteorological 
data  gathered  during  MIZEX-83;  to  report  findings  on  the 


degree  to  which  atmospheric 


anomalies 


that 


affect 


TABLE  1 


ELECTROMAGNETIC  SPECTRUM  SUMMARY 


TYPICAL 


FREQUENCY 

BAND 

FREQUENCY 

PRIMARY 

METHOD  OF 
PROPAGATION 

MAXIMUM 

DISTANCE  IN 
(NAUTICAL  MILES) 

Very  Low 
(VLF) 

3-30  KHz 

Ground  Wave 

5000 

Low  <LF> 

30-300  KHz 

Ground  Wave 

Sky  Wave 

1000-5000 

Medium 

(MF> 

300-3000  KHz 

Ground  Wave 

Sky  Wave 

100-1000 

1000-3000 

High  <HF> 

3-30  MHz 

Ground  Wave 

10-100 

100-250  (1  hop) 
100-1200  <>1  hop) 

Very  High 
(VHF) 

30-300  MHz 

Scatter 

Line  of  Sight 

800-1200 
lOO  (comm) 

300  (radar) 

Ultra-High 

(UHF) 

.3-3  GHz 

Scatter 

Line  of  Sight 

30-400 

100  (comm) 

300  (radar) 

Super-High 

<SHF> 

3-30  GHz 

Line  of  Sight 

100  (comm) 

300  (radar) 

Extremely- 
High  <EHF) 

30-300  GHz 

Line  of  Sight 

100 

Infrared 

CIR) 

1-400  THz 

Line  of  Sight 

5 

Optical 


400-800  THz  Line  of  Sight 


0-100 


subject  to  the  classical  dielectric  effects  described  in 
chapter  I I . 

The  troposphere  is  constantly  in  a  state  of  change  with 
different  layers,  air  masses,  wind  conditions,  and 
temperature  an  i  pressure  levels.  Generally,  humidity, 
temperature  and  pressure  change  with  altitude  and  so  does 
the  index  of  refraction.  It  will  be  shown  later  in  this 
chapter  how  the  index  of  refraction  can  be  redefined  in 
terms  of  humidity,  temperature  and  pressure.  Within  the 
troposphere,  layers  with  different  indices  of  refraction 
can  form  because  of  variations  in  humidity,  temperature  and 
pressure.  If,  because  of  changes  in  these  conditions, 
refraction  decreases  with  altitude  much  faster  than  normal, 
an  atmospheric  duct  can  be  formed.  Variations  in  humidity, 
temperature  and  pressure  sufficient  to  cause  a  duct  are  a 
direct  product  of  the  existing  weather  conditions. 

This  chapter  will  present  a  description  of  the 
atmospheric  parameters  and  mechanisms  involved  in  duct 
formation  and  will  describe  how  vertical  measurements  of 
tropospheric  humidity,  temperature  and  pressure  can  be  used 
to  calculate  refractive  gradients.  Finally,  it  will  be 
shown  how  analysis  of  these  gradients  can  provide 
predictions  of  atmospheric  effects  on  electromagnetic  wave 
propagation . 


A.  DUCTING  AND  REFRACTION 
1 .  Duct  Formation 

As  stated,  electromagnetic  waves  in  the  high 
frequency  band  can  be  refracted  back  to  the  Earth  by  the 
ionosphere.  These  waves  can  then  be  reflected  from  the 
Earth's  surface  back  to  the  ionosphere  to  be  refracted 
again.  As  this  process  is  repeated,  these  waves  can 
propagate  for  thousands  of  miles.  Electromagnetic  waves  in 
the  VHF,  UHF  and  SHF  frequency  bands  (see  table  1)  are 
unaffected  by  the  ionosphere.  However,  if  an 
electromagnetic  wave  in  one  of  these  frequency  bands 
encounters  a  gradient  of  refractive  index  and  is  within  a 
critical  angle,  it  can  be  refracted  back  toward  the  surface 
of  the  Earth.  In  this  way,  the  wave  can  be  channeled  along 
the  curvature  of  the  Earth.  This  region  is  called  an 
atmospheric  duct  and  the  preceding  case  was  an  example  of  an 
elevated  duct.  Duct  types  will  be  covered  in  detail  later 
in  this  chapter.  Electromagnetic  waves  in  a  duct  are 
directed  much  like  electromagnetic  waves  in  a  wave  guide  and 
ranges  can  be  extended  far  beyond  the  horizon  as  shown  in 
figure  6.  Ducts  are  primarily  caused  by  warm  dry  air 
overriding  cool  moist  air  resulting  in  a  temperature 
inversion.  It  can  be  seen  from  figure  7  (Helvey  and 
Rosenthal,  1983,  adapted  from  the  GTE  Sylvania  Report)  that 
this  happens  with  regularity  throughout  the  ocean  areas  of 
the  world. 


£7 


B 


FIGURE  6 

(a)  Radar  Wave  Propagation  under  "Standard'’ 
Atmospheric  Conditions,  (b)  Radar  Wave 
Propagation  Path  in  a  Surface  Duct. 
Reproduced  from  IREPS  Manual 
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Mean  Annual  Frequency  of  Ducting  <a)  Winter,  (b)  Summer 
From  GTE  Sylvania  Report 


2 


Ref ractivity 


The  index  of  refraction  can  be  scaled  to  yield 
atmospheric  refractivity  N.  N  is  related  to  n  by: 

N  =  <n-l>  x  106  (3.1) 

When  expressed  in  terms  of  the  atmospheric  parameters 
pressure  P  in  millibars,  temperature  T  in  degrees  Kelvin  and 
water  vapor  pressure  e  in  millibars,  N  becomes  (Kerr,  1951): 

N  =  t  <77 . S  P/T)  *  <3.73  x  105  e/T^)]  (3.2) 

Refractivity  can  be  directly  calculated  from 
humidity,  temperature  and  pressure.  But,  refractivity  at 
any  given  altitude  is  not  an  important  parameter  by  itself. 
A  vertical  profile  of  N  ,  N  values  plotted  against  altitude, 
or  a  plot  of  the  change  in  N  with  respect  to  the  change  in 
altitude  <dN/dz)  provides  a  graphic  method  by  which  a  duct 
can  be  identified.  However,  identifying  a  duct  and  defining 
its  boundaries  using  the  N  profile  requires  the  use  of  an 
overlay  and  is  a  cumbersome  technique.  To  make  the 
determination  of  duct  occurrence  more  convenient,  modified 
refractivity  has  been  developed. 

3.  Modified  Refractivity 

Modified  refractivity  M  takes  the  Earth's 
curvature  into  account  and  allows  quick  identification  of 


30 


ducting  whenever  M  values  decrease  with  altitude.  M  is 
related  to  N  by: 

M  =  N  ♦  < z/a )  X  106  (3.3) 

where  z  is  equal  to  altitude  in  kilometers  and  a  is  the  mean 


Earth 

radius 

in  kilometers.  Using  a  mean 

Earth 

radius 

of 

6.37 

X  103 

KM,  106/a  equals  157  KM-1. 

Now 

M  can 

be 

expressed  as 


M  =  N  ♦  157z  for  z  in  kilometers  (3.4) 

For  the  Standard  Atmosphere  where  there  are  no  well 
defined  layers  or  inversions,  temperature  and  pressure 
decrease  with  altitude.  As  shown  in  figure  8,  N  decreases 
and  M  increases  in  these  standard  conditions.  A  complete 
derivation  of  the  refractivity  and  modified  refractivity 
formulae  is  given  by  Kerr  (1951). 

4 .  Trapping 

A  trapping  layer  exists  when  dN/dz  is  less  than 
-157  or  dM/dz  is  less  than  zero.  Figure  9  shows 
representati ve  M  and  N  profiles  for  an  elevated  trapping 
layer.  Under  these  conditions  the  electromagnetic  wave 
will  be  bent  down  relative  to  the  Earth.  For  a  wave  to  be 
trapped,  a  ducting  condition  must  exist  where  the  wave  is 
trapped  between  the  upper  and  lower  boundaries  of  the  duct. 


Trapped  means  the  energy  transmitted  within  the  critical 


m  o  c  h 


FIGURE  8 


H  and  N  Profiles  for  Standard  Atmosphere 


acceptance  angle  will  be  partially  confined  within  the  duct 
with  very  little  loss.  In  the  atmosphere  n  varies  between 
1.0003  and  1.0005.  Using  equation  5.6  with  ni  equal  to 
1.0005  and  n2  equal  to  1.0003,  n2/ni  is  approximately  equal 
to  0.9998.  By  inserting  this  figure  in  equation  5.6,  the 
critical  angle  for  trapping  is  C90  -  ARCSINE  (0.9998)] 
degrees  which  is  approximately  equal  to  (90  -  89)  degrees 
or  one  degree.  Since  the  wave  front  is  assumed  to  be 
spherical,  the  entire  wave  front  will  not  meet  this 
requirement  and  not  all  of  the  energy  in  the  wave  will  be 
trapped.  Some  will  pass  through  the  layer  leading  to  the 
formation  of  radar  or  radio  holes  as  shown  in  figure  10. 
However,  for  the  wave  to  be  trapped,  the  emitter  must  be 
within  the  trapping  layer  to  meet  the  one  degree  critical 
acceptance  angle  requirement. 

Hitney  and  Paulus  (1979)  describe  refractive 
conditions  in  terms  of  the  refractive  gradient  found  in  the 
Standard  Atmosphere.  Under  standard  conditions,  values  for 
the  ft  gradient  range  between  79  and  157  ft  units/km.  Under 
these  conditions,  an  electromagnetic  wave  propagating  in  the 
troposphere  will  be  refracted  down,  but  with  a  curvature 
less  than  that  of  the  Earth.  Subrefraction  describes  an  M 
gradient  weaker  than  the  standard  gradient  or  greater  than 
157  ft  units/km.  Subrefracted  waves  will  be  refracted  upward 
away  from  the  Earth's  surface.  Superrefraction  describes  an 
II  gradient  stronger  than  the  standard  gradient  or  0  to  79  M 
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units/km.  Superref racted  waves  will  be  bent  downward  but  not 
trapped.  Trapping,  of  course,  occurs  when  the  M  gradient 
is  less  than  zero  or  M  decreases  with  height. 

B.  SURFACE  AND  ELEVATED  DUCTS 

Ducts  are  also  described  according  to  the  location  of 
their  lower  boundaries.  A  duct  with  its  lower  boundary  at 
the  surface  is  called  a  surface  duct;  however,  a  surface 
duct  may  extend  to  a  rather  high  altitude  <300  meters  to 
1  kilometer) .  A  duct  with  its  lower  boundary  above  the 
surface  is  called  an  elevated  duct.  The  thickness  of  the 
duct  is  determined  by  extending  a  vertical  line  down  from 
the  top  of  the  trapping  layer  (the  point  where  the  slope  of 
the  M  profile  changes  from  negative  to  positive)  to  the 
intersection  of  the  M  curve  as  shown  in  figure  11.  The  size 
of  the  duct  governs  frequencies  which  will  be  trapped. 
Basically,  the  thicker  the  duct,  the  lower  the  minimum 
trapped  frequency  will  be. 

C.  EVAPORATION  DUCT 

Evaporation  at  the  surface  of  the  sea  leads  to  near 
1005c  relative  humidity  in  the  air  immediately  above  the 
surface.  A  rapid  decrease  to  805c  to  905c  relative  humidity  a 
few  meters  above  the  surface  leads  to  a  refractive  gradient 
change  sufficient  to  bring  about  a  duct.  This  is  known  as 
the  evaporation  duct  and  exists  almost  all  the  time  over  all 
oceanic  areas. 
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As  shown  by  figure  12,  the  evaporation  duct  height, 
denoted  by  z«,  is  the  height  at  which  dM/dz  changes  from 
negative  to  positive  or  (dM/dz  =  0) . 

Paulus  and  Hitney  (1979)  have  found  that  typical  z* 
values  range  from  10  to  20  meters,  but  have  been  noted  as 
high  as  30  meters.  These  ducts  are  thinner  and  weaker  than 
other  ducts  and  trap  higher  frequencies.  However, 

significant  over-the-horizon  range  extensions  for  some 
radars  and  even  VHF/UHF  communications  can  be  realized. 

Evaporation  ducts  are  not  dependent  upon  the  large 
scale  meteorological  conditions  necessary  for  the  formation 
of  surface  and  elevated  ducts  and,  as  such,  will  not  be 
considered  in  this  analysis.  They  can  be  evaluated  from 
surface  temperature,  relative  humidity  and  wind  speed  data. 
Additional  information  on  evaluating  evaporation  ducts  is 
given  by  Paulus  and  Hitney  (1979). 

D.  MINIMUM  TRAPPED  FREQUENCY 

A  duct  of  thickness  “d”  can  be  compared  to  a  wave  guide 
of  a  specific  thickness.  The  wave  guide  is  manufactured  to 
exact  dimensions  according  to  the  proposed  operating 
frequency.  Likewise,  a  duct  will  trap  certain  frequencies 
depending  upon  its  thickness. 


FIGURE  12 

Evaporation  Duct  Refractive  Profile 
Adapted  from  IREPS  Manual 
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MIZEX-83 

DISTRIBUTION  BY  OPTIMUM  COUPLING  HEIGHT 


OPTIMUM  COUPLING  HEIGHT  (METERS) 


TABLE  2 


MIZEX-83  DUCT  ANALYSIS  DATA 


ASCENT 

NR 

OPTIMUM 

COUPLING 

HEIGHT 

DUCT 

THICKNESS 

(METERS) 

DUCT 

INTENSITY 
(M  UNITS) 

WIND  DIRECTION 
RELATIVE  ICE 
POSITION 

1 

8 

405.6 

111 

4.0 

260/0N 

ICE 

2 

10 

361.2 

77 

0.8 

180/0N 

ICE 

3 

13 

944.3 

216 

8.7 

190/0N 

ICE 

4 

14A 

392.7 

51 

1.2 

197/ON 

ICE 

5 

14B 

462.4 

84 

4.0 

197/ON 

ICE 

8 

15 

529.3 

34 

1.5 

224/ON 

ICE 

7 

16A 

795.8 

217 

17.3 

248/ON 

ICE 

8 

16B 

1717.1 

109 

6.0 

248/ON 

ICE 

9 

17A 

581.3 

112 

3.2 

212/ON 

ICE 

10 

17B 

2294.5 

191 

16.6 

212/ON 

ICE 

11 

18 

581.6 

113 

3.4 

227/ON 

ICE 

12 

20 

1488.5 

58 

O 

N 

158/ON 

ICE 

13 

22 

285.2 

81 

5.3 

077/OFF 

ICE 

14 

23 

307.8 

122 

10.4 

075/OFF 

ICE 

15 

34 

2096.8 

42 

1.9 

294/OFF 

ICE 

16 

67 

397.5 

31 

0.3 

300/0FF 

ICE 

17 

68 

2308.5 

69 

2.5 

230/0N 

ICE 

18 

77 

1820.4 

48 

0.4 

160/0N 

ICE 

19 

79 

1992.8 

53 

0.9 

155/ON 

ICE 

20 

85A 

575.1 

150 

12.5 

153/ON 

ICE 

21 

85B 

1219.9 

50 

0.8 

153/ON 

ICE 

22 

86 

1182.9 

123 

2.6 

100/0N 

ICE 

23 

88 

922.1 

76 

0.8 

160/0N 

ICE 

MEAN 

1028.84 

96.43 

4.58 

STND  DEV 

695.46 

54.73 

5.11 

5£ 


number,  for  each  sounding  where  a  duct  was  detected.  The 
first  plot  shows  the  overall  modified  refractivity  profile 
over  a  broad  altitude  range  and  the  second  concentrates  on 
the  area  of  the  duct.  In  four  of  the  soundings,  two  ducts 
at  different  altitudes  were  detected.  For  these  soundings, 
three  plots  were  made;  one  overall  profile  and  two  detailed 
profiles  with  "A"  indicating  the  lower  duct  and  "B”  indi¬ 
cating  the  upper  duct. 

C.  ANALYSIS 

Of  the  88  soundings  from  the  Polarbjorn  which  produced 
useful  meteorological  data,  19  contained  ducts.  Two  ducts 
were  present  in  four  of  these  soundings  bringing  the  total 
number  of  ducts  used  in  this  analysis  to  23.  Table  2  is  a 
list  of  optimum  coupling  height,  duct  thickness,  duct 
strength  and  the  mean  and  standard  deviation  for  each  para¬ 
meter  for  these  23  cases.  Table  2  also  lists  the  wind 
direction  in  degrees  and  assesses  the  wind  as  being  either 
"on  ice"  or  "off  ice"  for  each  of  the  cases.  This  is  a 
broad  assessment  based  on  the  Polarbjorn's  position  relative 
to  the  ice  pack  and  wind  direction.  These  parameters  will 
be  the  basis  for  establishing  a  refractivity  summary. 

1 .  Optimum  Coupling  Height 

Figure  17  is  a  frequency  distribution  by  optimum 
coupling  height  of  the  23  ducting  cases.  Figure  17  shows  14 


of  the  23  ducts  occurred  at  an  altitude  below  the  mean 


swing  in  the  PI  profile.  When  this  happens,  the  duct  may  be 


represented  as  stronger  than  it  actually  was.  In  the  worst 
case,  a  non-existent  duct  can  be  indicated  by  the  PI  profile. 
Some  MIZEX-83  radiosonde  data  is  questionable  for  this 
reason . 

B.  ANALYSIS  APPROACH 

A  mainframe  computer  was  used  for  initial  processing  of 
the  Polarbjorn's  radiosonde  sounding  data  tape.  Fortran 
routines  were  used  to  calculate  and  plot  refractivity  and 
modified  refractivity  profiles.  This  approach  was  adapted 
to  allow  for  relatively  quick  processing  of  the  voluminous 
amount  of  radiosonde  data.  Although  the  output  was  coarse, 
it  provided  a  means  for  initial  visual  identification  of 
ducts  using  the  modified  refractivity  profile.  Data  sets 
for  soundings  containing  ducts  were  then  processed  on  a 
desk-top  microcomputer  using  IREPS. 

Although  duct  height  and  thickness  can  be  read  directly 
from  the  IREPS  output,  a  more  detailed  plot  of  the  modified 
refractivity  profile  is  necessary  to  determine  duct  inten¬ 
sity.  To  obtain  the  necessary  plots,  modified  refractivity 
was  plotted  against  altitude  on  a  Zenith  Data  Systems  (Z- 
100)  desk-top  microcomputer.  Duct  optimum  coupling  height, 
duct  thickness  and  duct  intensity  were  read  directly  from 
these  plots  using  the  techniques  described  in  Chapter  III. 
Appendix  B  contains  two  plots. 


50 


with  corresponding  ascent 
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FIGURE  15 

Detailed  Course  Plot  for  the  Polarbjorn 
During  MIZEX-83 


Although  difficult  to  depict  on  figure  14,  the  Polarbjorn 
operated  within  or  in  close  proximity  to  these  ice  floes 
throughout  the  exercise.  Figure  15  is  a  more  detailed  plot 
of  the  Polarbjorn's  course.  These  two  figures  will  be 
referred  to  later  to  correlate  refractivity  conditions  to 
surface  meteorological  conditions  observed  during  the 
exercise . 

Of  the  103  radiosonde  ascents  made  from  the 
Polarbjorn,  88  provided  usable  data.  Appendix  A  lists  the 
time,  position  and  surface  observations  of  each  of  these  88 


soundings . 

F igure 

16 

is  a  position  plot 

of 

the 

88 

soundings 

and  can 

be 

related  to  figures  14 

and 

15 

for 

comparing  meteorological  and  refractivity  conditions. 

In  order  to  maximize  vertical  resolution,  the  ascent 
rates  were  kept  low  <120  to  150  meters/minute  versus  the 
normal  300  meters/minute.  With  temperature,  humidity  and 
pressure  being  measured  and  reported  every  10  seconds  and  a 
120  meter/minute  ascent, rate,  a  measurement  is  made  every  20 
meters . 

Although  radiosonde  soundings  generally  provide 
reliable  data,  some  conditions  can  lead  to  invalid  results. 
For  example,  surface  fog  can  saturate  the  humidity  sensor. 
In  this  condition,  high  humidity  level  readings  can  occur  in 
dry  air  well  above  the  fog.  When  this  sensor  abruptly  dries 
out,  actual  humidity  is  recorded  leading  to  an  abrupt  de¬ 
crease  in  the  humidity  profile  and  a  corresponding  negative 
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FIGURE  14 

Polarbjorn's  Course  and  Location  of  the  Ice  Pack 
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IV.  ACQUISITION.  ANALYSIS  AND  RESULTS 


A.  DATA  ACQUISITION 

1 .  Marginal  Ice  Zone  Experiment  <MIZEX> 

MIZEX  is  a  multi-year  meteorological  experiment 
designed  to  evaluate  meteorological  conditions  in  and  around 
the  Arctic  marginal  ice  zone  (MIZ) .  The  MIZ  is  a  zone  of 
roughly  100  km  width  characterized  by  the  transition  from 
packed  sea  ice  to  broken  ice  to  open  ocean.  MIZEX-83, 
conducted  during  July  and  August  of  1983  in  the  East  Green¬ 
land  Sea,  was  the  pilot  experiment  for  the  program.  Data 
were  gathered  using  meteorological  and  remote  sensing  air¬ 
craft  and  two  research  vessels  CPolarbjorn  and  Polarstern). 
Radiosondes,  sensors  attached  to  a  balloon  which  transmit 
measured  meteorological  data,  were  launched  from  the  re¬ 
search  vessels  at  various  points  in  the  MIZ.  The  radiosonde 
provides  the  information  necessary  to  plot  the  refractivity 
gradients  described  in  chapter  III.  This  analysis  will 


incorporate  only  radiosonde  data  taken  from  the  Polarbjorn, 


2 .  Radiosonde  Data 

Radiosonde  observations  on  the  Polarbjorn  were  made 
with  the  Vaisala  Micro-Cora  Upper  Air  sounding  system  using 
the  RS-80  radiosonde.  Figure  14  shows  the  course  of  the 
Polarbjorn  and  the  boundary  of  the  main  ice  pack.  In  sum¬ 
mer,  the  main  ice  pack  breaks  up  into  floes  at  the  edges. 


will  be  used  to  produce  an  area  refractivity  summary  based 
on  the  meteorological  data  gathered  during  MIZEX-83. 


» 


» 


provide  tactical  commanders  with  a  means  to  assess 
refractive  conditions  using  available  shipboard 

environmental  data.  Ultimately  IREPS  will  be  installed  on 
all  Navy  aircraft  carriers  and  used  as  an  aid  in  adjusting 
tactics  to  compensate  for  and  exploit  atmospheric  refractive 
effects,  (Beach,  1979) .  IREPS  processes  vertical 
temperature  and  humidity  profiles  on  the  Hewlett-Packard 
model  9845  desk-top  microcomputer  to  provide  four  basic 
product  outputs: 

1)  a  propagation  condition  summary  of  the  existing 
refractive  conditions  for  the  location  and  time 
of  the  input  data; 

2)  a  coverage  display  profile  of  a  specific  system 
based  on  the  systems  specifications  and  an  assess¬ 
ment  of  the  input  data; 

3)  a  computer  listing  of* the  input  data  set; 

4)  a  path  loss  display  representative  of  the 
performance  of  a  specified  system  for  a  given  data 
aet . 

If  on-scene  refractivity  data  are  not  available,  IREPS 
can  generate  products  that  show  probable  atmospheric  effects 
on  various  electronic  systems  using  a  historic  and 
geographical  refractivity  library.  Hitney  and  Paulus  (1979) 
give  a  detailed  description  of  IREPS  in  the  IREPS  Interim 
User's  Manual. 

IREPS,  designed  to  be  an  on-scene  tactical  refractivity 
assessment  aid,  was  used  to  analyze  part  of  the  data 
gathered  during  MIZEX-83.  In  the  following  chapter,  IREPS 
along  with  the  other  techniques  described  in  this  chapter 


marine  atmospheric  boundary  layer  <!1ABL>  .  Above  the  MABL, 
the  air  can  be  much  warmer  and  dryer.  At  the  boundary  of 
these  two  layers,  there  is  generally  a  sharp  decrease  in  N 
and  M.  Ducts  can  form  over  land  as  well  as  water  because 
the  variations  in  temperature,  pressure  and  humidity  with 
increasing  altitude  are  common  to  both  these  environments. 
Duct  persistence,  however,  is  enhanced  by  subsiding 
(sinking)  air  which  is  produced  by  trade  wind  circulation  of 
relatively  large,  stable  air  masses  in  the  area  20  to  30 
degrees  latitude  north  and  south  of  the  equator.  Trade  wind 
circulation  is  from  east  to  west  underneath  a  cap  of  warm 
dry  subsiding  air.  This  creates  a  trade  wind  inversion 
layer  and  the  positive  temperature  and  negative  humidity 
gradients  necessary  to  produce  ducting  conditions, 
(Neiburger,  et  al.,  1982).  These  weather  conditions 
occur  more  often  over  ocean  areas.  For  the  same  reason, 
the  persistence  of  ducting  conditions  in  the  extreme 
northern  and  southern  latitudes  is  rare.  Initial  duct 
formation  in  these  regions  is  based  on  smaller,  mesoscale 
meteorological  conditions  such  as  surface  flux  and  ice 
conditions . 

G.  ASSESSING  REFRACTIVE  CONDITIONS 

The  Naval  Ocean  Systems  Center  (NOSC) ,  San  Diego 
developed  and  continues  to  refine  the  Integrated  Refractive 
Effects  Prediction  System  (IREPS).  The  aim  of  IREPS  is  to 


REFRACTIVITY 


FIGURE  13 
Duct  Parameters 
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Based  on  experimental  values,  the  minimum  frequency  in 
Hertz  trapped  by  a  duct  of  thickness  d  is: 

fmin  =  3.6033  x  1011  x  cT<3/2)  C3.5) 

Kerr  <1951,  Chapter  2)  presents  a  detailed  examination  of 
duct  frequency  trapping. 

It  is  important  to  note  that  fmin  does  not  represent  a 
cut-off  frequency,  but  a  frequency  above  which  waves  will  be 
strongly  guided.  Frequencies  below  fmin  can  still  be 
affected  but  to  a  limited  degree. 

E.  DUCT  CHARACTERISTICS 

Knorr  <1982)  has  described  ducts  by  three  physical 
parameters- -duct  thickness,  intensity  and  optimum  couplinq 
height.  Duct  thickness  d  describes  the  vertical  size  of 
the  duct  and  relates  to  minimum  frequency  trapped.  Duct 
intensity  or  strength  is  described  in  M  units  and  is  the 
change  in  M  units  (A  ID  between  the  lower  boundary  of 
the  duct  and  optimum  coupling  height.  Optimum  coupling 
height  is  the  altitude  for  optimum  trapping  in  the  duct; 
altitude  where  the  critical  acceptance  angle  is  the  largest. 
Figure  13  shows  the  relationship  between  these  parameters. 

F.  METEOROLOGICAL  EFFECTS  LEADING  TO  DUCTING 

As  previously  stated,  duct  formation  is  dependent  on  the 
existence  of  two  differing  air  masses  or  layers.  Above  the 
ocean,  there  normally  exists  a  cool,  moist  layer  called  the 
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optimum  coupling  height  of  102S  meters  and  17  were  within 
one  standard  deviation. 

2.  Duct  Thickness 

Figure  18  is  a  frequency  distribution  by  duct  thick¬ 
ness  of  the  23  cases.  This  figure  shows  a  mean  thickness  of 
96.43  meters  with  19  of  the  23  cases  being  within  one  stan¬ 
dard  distribution. 

3 .  Duct  Intensity 

Figure  19  is  a  frequency  distribution  by  duct  inten¬ 
sity.  This  figure  shows  a  large  grouping  below  the  mean 
intensity  of  4.58  modified  refractivity  units.  Again,  19  of 
the  23  cases  are  within  one  standard  distribution. 

D.  RESULTS 

Duct  parameter  comparisons  made  in  this  section  are 
based  on  statistical  information  extracted  from  the  GTE 
Sylvania  report  (Ortenburgher ,  et  al.,  1978).  This  data 
base  contains  the  probability  of  occurrence  and  statistical 
parameters  such  as  thickness,  intensity  and  optimum  coupling 
height  of  elevated  ducts  collected  over  a  five  year  period 
at  sounding  stations  around  the  world.  For  this  analysis, 
comparisons  were  made  with  selected  mid  latitude  sounding 
stations . 

Figures  17,  18  and  19  reveal  the  following: 


(1)  Duct  optimum  coupling  height  in  the  MIZ  for  the 


MIZEX-83 

DISTRIBUTION  BY  DUCT  THICKNESS 


observation  period  does  not  vary  significantly  when  compared 
to  selective  mid  latitude  optimum  coupling  heights  of  1209 
meters  at  Birmingham,  Alabama,  1136  meters  at  Baghdad,  Iraq, 


and  727  meters  at  Qrendi,  Malta.  It  is  significant  to  note 
that  no  surface  based  ducts  were  detected. 

<2>  As  stated,  atmospheric  ducts  affect  electromag¬ 
netic  waves  in  the  UFH,  VHF  and  SHF  frequency  bands,  and 
duct  thickness  determines  minimum  trapped  frequency.  Using 
the  mean  duct  thickness  of  96  meters  and  applying  equation 
3.5  yields  a  minimum  trapped  frequency  of  approximately  383 
Mhz  which  is  in  the  UHF  band,  see  Table  1.  Applying  equa¬ 
tion  3.5  to  the  largest  and  smallest  detected  ducts  yields 
minimum  trapped  frequencies  of  approximately  110  Mhz  CVHF) 
and  2.1  Ghz  CSHF)  respectively.  However,  selected  mean  duct 
thicknesses  ranged  from  139  meters  in  Birmingham,  Alabama  to 
283  meters  in  Shreveport,  Louisiana;  all  significantly 
thicker  than  in  the  MIZ. 

<3)  Mean  duct  intensity  for  the  observation  period 
compares  favorably  with  many  mid  latitude  regions  of  the 
world.  However,  there  are  mid  latitude  regions  which 
average  much  stronger  ducts.  Compared  to  the  mean  intensity 
at  Shreveport,  Louisiana  of  7  M  units,  ducts  detected  during 
the  exercise  do  not  seem  very  weak.  However,  when  compared 
to  Kenitra  Morocco  with  a  mean  intensity  of  16  M  units. 


these  ducts  can  be  considered  weak  and  not  as  efficient  at 


trapping  and  channeling  energy  as  ducts  in  some  parts  of  the 
world . 

Inferences  about  the  existing  meteorological  conditions 
can  be  made  from  the  above  observations.  The  optimum  cou¬ 
pling  height  corresponds,  roughly,  to  the  height  of  the 
inversion  layer  and  weak  ducts  result  when  the  change  in 
the  refractive  index  is  small.  But,  to  get  a  change  in 
refractivity  requires  a  change  in  temperature  and/or  pres¬ 
sure  and/or  humidity.  A  major  question  to  be  answered  by 
this  analysis  is  what  observable  meterological  conditions  in 
the  MIZ  lead  to  ducting? 

Figure  20  is  a  histogram  of  duct  occurrence  in  time 
using  optimum  coupling  height  as  a  parameter.  The  time 
period  is  the  length  of  the  observation  period  or  the  dura¬ 
tion  of  MIZEX-83  scaled  to  the  launch  or  ascent  number. 
Figure  21  is  the  same  type  of  histogram  using  the  actual 
number  of  days  <36)  in  which  ascents  were  made  in  place  of 
the  ascent  number  and  reflects  the  time  period  more  accura¬ 
tely.  Since  these  histograms  are  designed  to  represent  only 
the  occurrence  or  non-occurrence  of  a  duct  at  some  point  in 
time,  only  one  of  the  ducts  in  each  of  the  four  cases  where 
double  ducts  occurred  has  been  considered.  The  determina¬ 
tion  of  which  duct  to  use  in  the  analysis  was  based  on  the 
height  of  the  stratus  cloud  layer  which  is  a  general  indica¬ 
tor  of  the  top  of  the  inversion  layer.  A  determination  as 
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to  what  meteorological  factors  caused  the  other  duct  in 
these  four  cases  has  not  been  made. 

Figures  20  and  21  clearly  show  that  ducts  occurred  in 
groups  and  that  there  were  periods  of  time  when  no  ducts 
were  detected.  Using  Table  2,  the  following  meteorological 
conditions  can  be  correlated  to  the  ducting  and  non-ducting 
periods  indicated  by  Figure  20  and  21. 

<1)  Ducting  conditions  were  prevalent  during  the 
period  1  to  4  July.  Eleven  soundings  during  these  four  days 
detected  nine  ducts.  Ducts  were  also  previously  detected  on 
27  and  28  June.  During  this  period,  the  Polarbjorn  was  in 
broken  ice  floes  at  the  edge  of  the  pack.  Winds  were  light 
and  predominantly  from  the  south  south-west  onto  the  ice. 
Fog  and  stratus  overcast  dominated  the  period.  Referring  to 
Figure  16,  these  ducts  were  detected  by  the  group  of  sound¬ 
ings  indicated  in  the  upper  right  corner  of  the  radiosonde 
ascent  plot. 

<2>  From  5  to  19  July,  43  soundings  detected  only 
one  duct.  As  the  Polarbjorn  moved  south  and  away  from  the 
main  ice  pack,  ducts  were  no  longer  present.  Meteorological 
conditions  during  this  period  were  also  characterized  by 
stratus  overcast  and  fog.  Winds  remained  light  and  after  7 
July,  blew  predominantly  from  the  north  off  the  ice.  The 
soundings  taken  during  this  period  are  the  group  in  the 
right  center  portion  of  Figure  16. 


(3)  Figure  20  shows  another  group  of  ducts  that 
were  detected  during  the  period  20  to  29  July  with  the 
majority  occurring  after  23  July.  From  23  to  29  July  twelve 
soundings  detected  five  ducts.  The  Polarbjorn  was  east  of 
broken  ice  floes  during  this  period  with  the  wind  blowing 
from  the  east  and  south  onto  the  ice.  This  period  was  also 
characterized  by  fog  and  stratus  overcast.  This  group  of 
soundings  can  be  seen  in  Figure  16  in  the  lower  left  portion 
of  the  plot. 


SUMMARY  AND  RECOMMENDATIONS 


A.  SUMMARY 

The  optimum  coupling  height  and  thickness  of  ducts  de¬ 
tected  during  MIZEX-83  do  not  vary  significantly  from  ducts 
found  in  other  parts  of  the  world.  However,  these  ducts  do 
not  tend  to  be  as  intense  as  ducts  found  in  some  mid  lati¬ 
tude  locations.  In  some  cases,  radiosonde  data  collected 
may  not  accurately  reflect  refractive  conditions  because  of 
a  fog  saturated  humidity  sensor.  This  does  not  necessarily 
mean  that  the  ducts  did  not  exist,  but  that  they  may  be 
represented  as  stronger  than  they  actually  were.  L.N. 
Orthenburger  (NOSC  TD  260,  1979)  states  that  radiosonde  data 
tend  to  underestimate  ducting  occurrence  because  of  time- 
lags  and  other  inaccuracies  in  the  humidity  sensing  element. 
This  characteristic  of  the  radiosonde  may  partially  offset 
the  inaccuracies  caused  by  saturation,  but  may  cause  the 
rest  of  the  ducts  to  be  represented  as  weaker  than  they 
actually  were. 

During  the  period  of  this  experiment  ducting  conditions 
occurred  significantly  more  often  when  the  wind  blew  from 
the  open  sea  onto  the  ice.  Based  on  this  study,  no  other 
general  or  specific  meteorological  conditions  have  been 
correlated  to  the  ducting  and  non-ducting  periods.  Light 
winds,  fog,  stratus  overcast,  drizzle  and  proximity  to  the 
ice  were  characteristic  of  the  experiment.  Although  no 
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ducts  were  detected  during  the  brief  clear  periods,  the  data 


set  is  too  small  to  directly  correlate  these  two  conditions. 

B.  ELECTRONIC  WARFARE  ASSESSMENT 

Of  importance  to  a  tactical  commander  operating  in  the 
MI2  would  be  the  ability  to  assess  atmospheric  refractivity 
conditions  to  gain  a  tactical  electronic  warfare  advantage. 
Of  course,  this  assessment  is  made  possible  by  real  time 
assessment  aids  such  as  IREPS,  but  a  knowledge  of  typical 
area  conditions  would  enhance  mission  planning.  For  exam¬ 
ple,  knowledge  of  the  large  scale  ducting  conditions  in  the 
Indian  Ocean  allows  a  tactical  commander  to  stand  off  and 
monitor  adversary  VHF/UHF  communications  at  a  range  well 
over  the  horizon.  Since  this  analysis  is  a  first  step  in 
analyzing  and  documenting  refractivity  conditions  in  the 
niZ,  only  general  statements  about  electronic  warfare  opera¬ 
tions  con  be  mode.  Therefore,  the  following  observations 
should  be  considered  in  the  context  of  the  size  of  the  data 
set  and  the  time  of  year  in  which' the  data  were  taken. 

The  lack  of  surface  based  ducts  will  restrict  over-the- 
horizon  (OTH)  electronic  support  measures  (ESM)  in  the 
VHF/UHF/5HF  frequency  range  from  surface  platforms. 
However,  some  OTH  ESM  may  be  possible  if  an  evaporation  duct 
is  present.  Likewise,  if  this  is  the  case,  a  task  force 
could  propagate  electromagnetic  energy  in  these  frequency 
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ranges  without  fear  of  OTH  detection  by  an  adversary 
operating  on  the  surface. 

Airborne  ESP!  as  well  as  electronic  countermeasures  (ECM) 
such  as  jamming,  range  gate  pull-off  and  false  target  gen¬ 
eration  could  be  enhanced  because  of  the  existence  of 
elevated  ducts.  Extended  ranges  could  be  exploited  as  well 
as  radar  holes  by  knowledge  of  the  existence  of  and  optimum 
coupling  height  of  MIZ  ducts.  These  ducts  do  not  tend  to  be 
strong  and  seem  dependent  on  an  on-ice  wind.  Further,  if 
this  wind  is  caused  by  the  sea-ice  transition,  ducts  may  not 
extend  far  from  the  boundary  of  the  ice  pack.  These  factors 
limit  the  possible  tactical  advantage.  As  in  the  case  of 
surface  operations,  the  opposite  also  exists  for  airborne 
platforms;  PIIZ  ducts  are  capable  of  providing  extended 
detection  and  exploitation  ranges  for  an  adversary. 

C.  RECOMMENDATIONS 

This  analysis  has  been  limited  by  the  size  of  the  data 
set.  The  following  recommendations  address  areas  of 
analysis  that  were  not  considered  or  covered  in  sufficient 
detail . 

1)  The  mesoscale  meteorological  conditions  leading  to 
ducting  in  the  MIZ  were  only  briefly  mentioned  here. 
However,  this  is  the  thesis  of  another  analysis  effort  that 
will  correlate  MIZEX-83  radiosonde  data  with  other  data  such 
as  satellite  photographs.  This  analysis  should  provide  more 


insight  into  the  cause  of  the  second  duct  in  the  double 
ducting  cases  and  the  actual  effect  of  surface  wind  on  PIIZ 
duct  formation. 

2)  Sufficient  data  were  not  available  to  accurately 
determine  the  spatial  extent  of  the  ducts  detected  during 
MIZEX-83.  All  soundings  were  taken  in  the  proximity  of  the 
ice,  therefore,  it  was  not  possible  to  determine  if  the 
ducts  extended  over  the  ice,  the  open  water  or  both. 
Further,  it  could  not  be  determined  how  far  the  ducts 
extended  from  the  ice  edge.  Knowledge  of  the  spatial  extent 
of  MIZ  ducts  would  be  extremely  important  in  any  electronic 
warfare  planning.  Data  gathered  in  future  MIZ  experiments 
should  take  this  into  account  and  take  soundings  in  open 
water  at  a  distance  from  the  ice. 

3)  Knorr  (1982)  proposes  that  there  is  a  correlation 
between  duct  thickness  and  intensity  and  has  presented 
methods  to  compute  the  power  distribution  of  signals 
propagated  through  elevated  ducts  using  the  joint 
probability  of  elevated  duct  parameters.  An  analysis  of 
this  type  on  MIZEX-83  data  would  be  valuable  in  determining 
expected  range  improvements  when  propagating  in  a  duct. 

4)  MIZEX-84  has  just  recently  been  completed  providing 
another  data  set  which  will  reveal  more  information  about 
the  area.  This  should  be  added  to  MIZEX-83  results  to 
increase  the  knowledge  of  refractive  conditions  in  the  MIZ. 
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MIZEX-83  RADIOSONDE  ASCENTS  DATA 


ASNT  TIME  TEMP  WIND 

NR  DATE  LAT  LONG  GMT  PRES  CNTG  DIR 


1 

6/23 

80.08N 

4.02E 

2330 

999 . 1 

3.1 

180 

STRATUS,  IX 

ICE,  WATER  l.OC 

2 

6/24 

80.60N 

4. 10E 

1136 

996 . 8 

0.9 

148 

STRATUS,  LT 

SNOW,  IN 

1  ICE, 

CEILING 

200 

METERS 

3 

6/24 

80.86N 

5.37E 

2339 

1003.5 

0.7 

194 

STRATUS,  FOG, 

IN  ICE 

4 

6/25 

80.87N 

5.58E 

1139 

1006.2 

1.2 

120 

STRATUS,  NO  WIND 

5 

6/25 

81 .OON 

5.30E 

2333 

1001.0 

1.5 

154 

FOG, 

DRIZZLE 

6 

6/26 

81 .02N 

5.56E 

1135 

999.6 

0.0 

210 

THIN 

FOG,  IN 

ICE 

7 

6/27 

81 .07N 

6. 12E 

1138 

1000.1 

-0.7 

220 

THIN 

FOG,  BLUE  SKY  OVERHEAD 

.  IN  ICE 

8 

6/27 

81 .05N 

6.07E 

2331 

1000.2 

-1.0 

260 

FOG 

9 

6/28 

81 .05N 

6.07E 

1135 

1002.3 

-2.1 

330 

LOW  THIN  OVERCAST 

10 

6/28 

81 .07N 

6.  HE 

2331 

1011.9 

-3.7 

180 

NONE 

11 

6/29 

81 . 12N 

6.30E 

2340 

1012.1 

-0.8 

194 

CLEAR 

12 

6/30 

81 . 17N 

6.86E 

1133 

1012.8 

0.7 

235 

STRATUS,  CEILING  100 

METERS 

13 

7/01 

81.17N 

6.80E 

0023 

1013.2 

0.5 

190 

STRATUS 

14 

7/01 

81 .20N 

7.07E 

0943 

1010.1 

-0.7 

197 

THIN 

FOG,  BLUE  SKY  OVERHEAD 

83 

7/25  79.67N  2.98E 

2332 

1009.0 

1.6 

200 

6.0 

STRATOCUMULUS 

84 

7/26  79.02N  2.65E 

1937 

1016.4 

0.1 

240 

3.0 

STRATUS,  CEILING  80 

METERS 

;v  \- 

85 

7/27  78.57N  1.96E 

1404 

1019.5 

4.1 

153 

5.4 

*  • 

■*.  *"•  V 

STRATOCUMULUS,  IN  OPEN  WATER,  ICE 

TO  THE 

WEST 

86 

7/28  77.99N  -3.45E 

1331 

1013.1 

1.8 

100 

O 

in 

STRATOCUMULUS  WITH 

BREAKS 

AND  OCCASIONAL 

DRIZZLE 

87 

7/28  77.99N  -4.32E 

1737 

1010.1 

1.8 

102 

7.6 

HEAVY  FOG 


88  7/29  77.18N  -5.87E  1329  1003.5  0.6 

STRATUS,  CEILING  400  METERS 


160 


4.0 


7/19  79.17N  0.88E  2327  1010.3  -0.1 

STRATUS,  LT  SNOW  OR  SLEET 


340 


7/20 

FOG 

79.00N 

-1 .07E 

1127 

1013.5 

-0.3 

300 

7/20 

78.97N 

-1 .39E 

1805 

1013.1 

0.0 

230 

MOSTLY  CLEAR, 

1/8  ALTOCUMULUS  TO  THE  WEST 

7/20 

78.97N 

-1 .56E 

2329 

1013.1 

-0.3 

230 

SUNNY 

,  CLEAR  • 

7/21 

CLEAR 

78.90N 

-1 .70E 

0450 

1013.6 

-0.3 

185 

7/21 

78 . 98 N 

-1.90E 

1131 

1014.8 

1 . 1 

200 

SUNNY 

,  ALTOCUMULUS  AND  ALTOSTRATUS 

TO  THE 

SOUTH 
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